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Abstract 

The sequence-dependent folding landscapes of nucleic acid hairpins reflect much of the complexity of 
biomolecular folding. Folding trajectories, generated using single molecule force clamp experiments by 
attaching semiflexible polymers to the ends of hairpins have been used to infer their folding landscapes. 
Using simulations and theory, we study the effect of the dynamics of the attached handles on the handle- 
free RNA free energy profile F°Jz m ), where z m is the molecular extension of the hairpin. Accurate 
measurements of F° {z m ) requires stiff polymers with small L/l p , where L is the contour length of the 
handle, and l p is the persistence length. Paradoxically, reliable estimates of the hopping rates can only 
be made using flexible handles. Nevertheless, we show that the equilibrium free energy profile F° q (z m ) 
at an external tension f m , the force (/) at which the folded and unfolded states are equally populated, 
in conjunction with Kramers' theory, can provide accurate estimates of the force-dependent hopping 
rates in the absence of handles at arbitrary values of /. Our theoretical framework shows that z m is a 
good reaction coordinate for nucleic acid hairpins under tension. 
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A molecular under standingof how proteins and RNA fold is needed to describe the functions 



s [l| a 



2] , interactions between biomolecules, and the origins of misfolding 



that is linked to a number of diseases jjj. The energy landscape perspective has provided a 
conceptual framework for describing the mechanisms by which unfolded molecules navigate the 



large conformational space in search of the native state 



a a. 



Recently, single molecule 



techniques have been used to probe features of the energy landscape of proteins and RNA that 



are not easily accessible in ensemble experiments 
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It is possible to construct the shape of the energy landscape, including the energy scales of 
ruggedness Q, Q, using dynamical t—ej .that are generated by applying a constant 
force (/) to the ends of proteins and RNA [5|, [14J, [21|, |22|. If the observation time is long 



enough for the molecule to sample the accessible conformational space, then the time average of 
an observable X recorded for the a th molecule ((X) = lim^oo | f Q dTX a {T)) should equal the 
ensemble average {{X) = lim^^oo j^J2f =1 Xi), and the distribution P(X) should converge to 
the equilibrium distribution function P eq (X). Using this strategy, laser optical tweezer (LOT) 
experiments have been used to obtain the sequence-dependent folding landscape of a number of 
RNA and DNA hairpins (s, 2, Il4 , |23| . using X = R m , the end-to-end distance of the hairpin 



that is conjugate to /, as a natural reaction coordinate. In LOT experiments, the hairpin is 



]), whose ends are attached 



held between two long handles (DNA |5|] or DNA/RNA hybrids 
to polystyrene beads (Fig. la). The equilibrium free energy profile (3F eq (R m ) = — log P eg (R m ) 
{(3 = 1/ksT, ks is the Boltzmann constant, and T is the absolute temperature) may be useful 
in describing the dynamics of the molecule, provided R m is an appropriate reaction coordinate. 

The dynamics of the RNA extension in the presence of / (z m = Zy — ~ R m , provided 
transverse fluctuations are small) is indirectly obtained in an LOT experiment by monitoring 



the distance between the attached polystyrene beads [z, 



sys 



one of which is optically 



trapped at the center of the laser focus (Fig la). The goal of these experiments is to extract the 
folding landscape (f3F°Jz m )) and the dynamics of the hairpin in the absence of handles, using 
the /-dependent trajectories z sys (t). To achieve these goals, the fluctuations in the handles 
should minimally perturb the dynamics of the hairpin in order to probe the true dynamics of 
a molecule of interest. However, depending on L and l v (L is the contour length of the handle 
and l p is its persistence length), the intrinsic fluctuations of the handles can not only disort 
the signal from the hairpin, but also directly affect its dynamics. The first is a problem that 
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pertains to the measurement process, while the second is a problem of the coupling between the 
instruments and the dynamics of RNA. 

Here, we use coarse-grained molecular simulations of RNA hairpin and theory to show that, 
in order to obtain accurate /3F eq (R m ), the linkers used in the LOT have to be stiff, i.e., the 
value of L/lp has to be small. To investigate the handle effects on the energy landscape and 
hopping kinetics, we simulated the hairpin dynamics under force-clamp conditions by explicitly 
modeling the linkers as polymers with varying L and l p . Surprisingly, the force-dependent 
folding and unfolding rates that are directly measured using the time traces, z m (t), are close 
to the ideal values (those that are obtained by directly applying /, without the handles, to 
the 3' end with a fixed 5' end) only when the handles are flexible. Most importantly, accurate 
estimates of the /-dependent hopping rates over a wide range of /-values, in the absence of 
handles, can be made using j3F eq (R), in the presence of handles obtained at / = f m , the 
transition midpoint at which the native basin of attraction (NBA) and the unfolded basin of 
attraction (UBA) of the RNA are equally populated. The physics of a hairpin attached to 
handles is captured using a generalized Rouse model (GRM), in which there is a favorable 
interaction between the two non-covalently linked ends. The GRM gives quantitative agreement 
with the simulation results. The key results announced here provide a framework for using 
the measured folding landscape of nucleic acid hairpins at / ~ f m to obtain f -dependent folc - 



ing and unfolding times and the transition state movements as / is varied 
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RESULTS and DISCUSSIONS 



Modeling the LOT experiments: In order to extract the folding landscape from LOT 



experiments, the time sea 
have to be well-separated 



es associated with the dynamics of the beads, handles, and the hairpin 
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331 ] . The bead fluctuations are described by the overdamped 



Langevin equation ^dx p /dt = —kx p + F(t) where k is the spring constant associated with the 
restoring force, and the random white-noise force F(t) satisfies (F(t)) = and (F(t)F(t')) = 
2 r jkBT5(t — t'). The bead relaxes to its equilibrium position on a time scale r r = j/k. In terms of 
the trap stiffness, k p , and the stiffness k m associated with the Handle-RNA-Handle (H-RNA-H; 
see Fig. 1), k — k p + k m . With 7 = 6nria, a = 1/xm, 77 
k m ~ 0.1 pN/nm, we find r r < 1 ms. In LOT experiments 
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cP, k„ « 0.01 pN/nm |34j, and 
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33| , separation in time scales 
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is satisfied such that rfc ~ tJ. 3> r r at / ~ f m , where and rf. are the intrinsic values of the 
RNA (un) folding times in the absence of handles. 

Since Zfy^ IS 3i natural reaction coordinate in force experiments, the dispersion of the bead 
position may affect the measurement of F eq (z m ). At equilibrium, the fluctuations in the bead 
positions satisfy 5x 2 eq ~ ksT/ikp + k m ) ~ ksT/km, and hence k m should be large enough to 
minimize the dispersion of the bead position. The force fluctuation, 5f 2 q ~ kBTk 2 /{k p + k m ), 
is negligible in the LOT because k p <C k m , and as a result Sf eg / f m ~ 0, since 5f eq ~ 0.1 pN 
while / m ~ 15 pN. Thus, we model the LOT setup by assuming that the force and position 
fluctuations due to the bead are small, and exclusively focus on the effect of handle dynamics 
on the folding landscape and hopping kinetics of RNA (Figs. la-b). 



Short, stiff handles are required for accurate free energy profiles: For purposes of 
illustration, we used the self-organized polymer (SOP) model of the P5GA hairpin 35j, and 
applied a force / = f m ~ 15.4 pN. The force is exerted on the end of the handle attached to 
the 3' end of the RNA (P in Fig. la), while fixing the other end (O in Fig. la). Simulations of 
P5GA with handles of length L = 25 nm and persistence length l p = 70 nm show that the 
extension of the entire system (z sys = z p — zq) fluctuates between two limits centered around 
z sys m 50 nm and z sys « 56 nm (Fig. lb). The time-dependent transitions in z sys between 50 
nm and 56 nm correspond to the hopping of the RNA between the Native Basin of Attraction 
(NBA) and Unfolded Basin of Attraction (UBA). Decomposition of z sys as z sys = z^ + z m + z^, 
where z%(= zy — z Q ) and z%(= z p — zy) are the extensions of the handles parallel to the force 
direction (Fig. la), shows that z sys (t) reflects the transitions in z m (t) (Fig. lb). Because the 
simulation time is long enough for the harpin to ergodically explore the conformations between 
the NBA and UBA, the histograms collected from the time traces amount to the equilibrium 
distributions P eq {X) where X = z sys , z%, z m , or z^ (Fig lb; for P eq (z^) and P eq (z^), see the 
Supporting Information SI Fig. 6a). To establish that the time traces are ergodic, we show that 
z T (t) = | Jq dTz sys (r) reaches the thermodynamic average (~ z sys P eq (z sys )dz sys =53.7 nm) 
after t > 0.1 sec (the magenta line on z sys (t) in Fig. lb). 

Fig. lb shows that the positions of the handles along the / direction fluctuate, even in 
the presence of tension, which results in slight differences between P eq (z sys ) and P eq (z m ). 
Comparison between the free energy profiles obtained from the z sys (t) and z m (t) can be used to 



1 



investigate the effect of the characteristics of the handles on the free energy landscape. To this 
end, we repeated the force-clamp simulations by varying the contour length (L = 5 — 100 nm) 
and persistence length (l p = 0.6 and 70 nm) of the handles. Fig. 2 shows that the discrepancy 
between the measured free energy F eq (z sys ) (dashed lines in blue) and the molecular free energy 
F e q(z m ) (solid lines in red) increases for the more flexible and longer handles (see the SI text 
and SI Fig. 6 for further discussion of the dependence of the handle fluctuations on L and l p ). 
For small l p and large L, the basins of attraction in F eq (z m ) are not well resolved. The largest 
deviation between F eq (z sys ) and F eq (z m ) is found when l p = 0.6 nm and L = 25 nm (L/l p m 40) 
(the graph enclosed by the orange box in Fig. 2a). In contrast, the best agreement between 
F e q(z S ys) and F eq (z m ) is found for l p = 70 nm and L = 5 nm (the graph inside the magenta 
box in Fig. 2), which corresponds to L/l p ~ 0.07. In the LOT experiments, L/l p ~ 6—7 [ol B. \h\ . 

Generalized Rouse model (GRM) captures the physics of H-RNA-H under ten- 
sion: In order to establish the generality of the relationship between the free energy profiles 
as measured by z m and those measured by z sys , we introduce an exactly solvable model that 
minimally represents the RNA and handles (Fig. 3a). We mimic the hairpin using a Gaussian 
chain with Nq monomers and Kuhn length a. The endpoints of the RNA mimic are harmonically 
trapped in a potential with stiffness k as long as they are within a cutoff distance c = 4nm. Two 
handles, each with monomers and Kuhn length b, are attached to the ends of the RNA (see 
Methods). We fix one endpoint of the entire chain at the origin, and apply a force f m ~ 15.4 pN 
to the other end. The free energies as a function of both the RNA's extension, R m = |r 3 / — r 5 /| 
(^ z m at high /) and the system's extension R sys = \rp — r | (~ z sys at high /) are exactly 
solvable in the continuum representation. We choose k such that f m is near the midpoint of the 
transition, so that d 3 r P eq (r) ps J°° d 2 r P eq (r). We tune iVo so that the barrier heights for the 
GRM and P5GA are similar at / = f m . These requirements give N = 20 and k m 0.54 pN/nm. 

While the stiffness in the handles of the simulated system (Fig. 1) cannot be accurately 
modeled using a Gaussian chain, the primary effect of attaching the handles is to alter the fluc- 
tuations of the endpoints of the RNA. By equating the longitudinal fluctuations for the WLC, 
(5H^)wlc ~ Lip 1/,2 (/5/) -3 ^ 2 , with the fluctuations for the Gaussian handles, (5~R^)g ~ Lb, we 
estimate that the effective persistence length of the handles scales as Ifff ~ b~ 2 f~ 3 (see the 
SI for details). Thus, smaller spacing in the Gaussian handles in the GRM will mimic stiffer 
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handles in the H-RNA-H system. The free energies computed for the GRM, shown in Fig. 3b-c, 
are consistent with the results of the simulations. The free energy profiles deviate significantly 
from F° (z m ) as increases or 'stiffness' decreases. The relevant variable that determines 
the accuracy of F eq (z sys ) is N^b 2 ~ Ljl e J^ with the free energies remaining unchanged if N^b 2 
is kept constant. The barrier height and well depths as a function of z m are unchanged as a 
function of L and b. However, the apparent activation energy is decreased as measured by 
z sys (seen in Fig. 2 as well). The GRM confirms that accurate measurement of the folding 
landscape using z sys requires stiff handles. 

Accurate estimates of the hopping kinetics requires short and flexible handles: 

Because LOT experiments can also be used to measure the force-dependent rates of hopping 
between the NBA and the UBA, it is important to assess the influence of the dynamics of the 
handles on the intrinsic hopping kinetics of the RNA hairpin. In other words, how should the 
structural characteristics of the linkers be chosen so that the measured hopping rates using the 
time traces z(t) and the intrinsic rates close as possible? 

Folding and unfolding rates of P5GA and the free energy profile without handles : We first 
performed force clamp simulations of P5GA in the absence of handles to obtain the intrinsic (or 
ideal) folding (t f ) or unfolding {jy) times, that serve as a reference for the H-RNA-H system. 
To obtain the boundary conditions for calculating the mean refolding and unfolding times, we 
collected the histograms of the time traces and determined the positions of the minima of the 
NBA and UBA, zp = 1.9 nm and Zy = 7.4 nm (Fig. 4a). The analysis of the time traces 
provides the transition times in which z m reaches z m = zy starting from z m = zp. The mean 
unfolding time Ty is obtained using either Ty = ^-/N^2 i Ty(i), or from the fits to the survival 
probability P F (t) = e"'/ 7 " 17 (SI Fig. 8). The mean folding time is similarly calculated, and the 
two methods give similar results. The values of t v and r F computed from the time trace of 
z m (t) are 2.9 ms and 1.9 ms, respectively. At f m = 15.4 pN and L = nm, the equilibrium 
constant K eq = t f /t v = 0.67, which shows that the bare molecular free energy is slightly tilted 
towards NBA at / = 15.4 pN. 

Hopping times depend on the handle characteristics : The values of the folding (rj? 1 ) and 
unfolding (r™) times were also calculated for the P5GA hairpin with attached handles (Fig. 1). 
As the length of the handles increases both Ty and t™ increase gradually, and the equilibrium 
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distribution shifts towards the UBA, i.e. K eq = t™ jr™ increases (Fig. 4b). Strikingly, the use of 
flexible handles results in minimal deviations of t™ and T™ from their intrinsic values (Fig.4b). 
Attachment of handles (stiff or flexible) to the 5' and 3' ends restricts their movement, which 
results in a decrease in the number of paths to the NBA and UBA. Thus, both t™ and 
increase (Fig. 4b). As the stiffness of the handle increases the extent of pinning increses. These 
arguments show that flexible and short handles, that have the least restriction on the fluctuations 
of the 5' and 3' ends of the hairpin, cause minimal perturbation to the intrinsic RNA dynamics, 
and hence the hopping rates. 

Because the experimentally accessible quantity is the extension of the H-RNA-H, it is impor- 
tant to show that the transition times can be reliably obtained using z sys (t). Although z sys (t) 
differs from z m (t) in amplitude, the "phase" between the two quantities track each other reli- 
ably throughout the simulation, even when the handles are long and flexible (see SI Fig. 7). We 
calculated t^ s and T s F ys by analyzing the trajectories z sys (t) using the same procedure used to 
compute their intrinsic values. Comparison of Ty S (t f vs ) and t™ (t™) for both stiff and flexible 
handles shows excellent agreement at all L values (Fig.4b). Thus, it is possible to infer the RNA 
dynamics z m {t) by measuring z sys (t). 

Theoretical predictions using the GRM are consistent with the simulations: The simulation 
results can be fully understood using the GRM (Fig 3a), for which we can exactly solve the 
overdamped Langevin equation using the discrete representation of the Gaussian chain (see 
Methods). By assuming that transverse fluctuations are small (which is reasonable under the 
relatively high tension of / = 15.4 pN), we use the Wilemski and Fixman (WF) theory 36] 
to determine an approximate time of contact formation (r™ = (&™) -1 ) as a function of b (i.e. 
increasing handle 'stiffness') and A^. The refolding rate of the RNA hairpin under tension is 
analogous to A;™. A plot of /c^(L)//c™(0) versus L (Fig. 4c) illustrates that smaller deviations 
from the handle-free values occur when l p is small. Moreover, Fig. 4c shows that the refolding 
rate decreases for increasing regardless of the stiffness of the chain. The saturating value of 
k™ as Nh — > oo depends on b, with 'stiffer' handles having a much larger effect on the folding 
rate. While the handles used in LOT experiments are significantly longer than the handle lengths 
considered here, the saturation of the folding rate suggests that L ~ 100 nm is sufficiently long 
for finite-size effects to be negligible. 

We also find the dependence of hp on L agrees well with the behavior observed in the 
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simulation of P5GA. The ratio k'p(L) /fc^(O) for b = a agrees well with the trends of the flexible 
linker (l p = 0.6 nm) for all of the simulated lengths, with both saturating at kp{L) ps 0.35^(0) 
for large L. The trends for 'stiffer' chains (smaller b) in the GRM qualitatively agree with 
the P5GA simulation with stiff handles (l p = 70 nm), with remarkably good agreement for 
0.1 < b/a < 0.2 over the entire range of L. The GRM, which captures the physics of both the 
equilibrium and kinetic properties of the more complicated H-RNA-H, provides a theoretical 
basis for extracting kinetic information from experimentally (or computationally) determined 
folding landscapes. 



Free energy landscapes and hopping rates: Stiff handles are needed to obtain F eq (z sys ) 
that resembles F° q (z m ), whereas the flexible handles produce hopping rates that are close 
to their handle-free values. These two findings appear to demand two mutually exclusive re- 
quirements in the choice of the handles in LOT experiments. However, if z m is a good reaction 
coordinate, then it should be possible to extract the hopping rates using accurately measured 
F eq (z sys )(tt F eg (z m ) ps F° q (z m )) at / ~ f m , using handles with small L/l p . The (un)folding times 
can be calculated using the mean first passage time (Kramers' rate expression) with appropriate 



boundary conditions 37] 
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where z min , z max , z v and zp are defined in Fig. 4a. The effective diffusion coefficient D f (Djj) 
is obtained by equating Tp R {jy R ) in equation (jTJ), with F eq (z m ) = F° q (z m ), to the simulated 
t f ( T u)- We calculated the /-dependent T™(f) and r™(f) by evaluating equation (jTJ) using 
F° q (z m \f) = F° q (z m \f m ) — (/ — fm) ■ z m . The calculated and simulated results for P5GA are 
in good agreement (Fig 5a-b). At the higher force (/ = 16.8 pN), the statistics of hopping 
transition within our simulation time is not sufficient to establish ergodicity. As a result, the 
simulation results are not as accurate at high forces (see SI Fig. 9). To further show that the use 
of F° Jz m \f) in equation (JTJ) gives accurate hopping rates, we calculated r^(/) for the GRM and 
compared the results with direct simulations of the handle-free GRM, which allows the study of 
a wider range of forces (see Methods). The results in Fig. 5c show that F°Jz m \f) indeed gives 
very accurate values for the transition times from the UBA and NBA over a wide force range. 

8 



CONCLUSIONS 



The self-assembly of RNA and proteins may be viewed as a diffusive process in a multi- 
dimensional folding landscape. To translate this physical picture into a predictive tool, it is 
important to discern a suitable low-dimensional representation of the complex energy landscape, 
from which the folding kinetics can be extracted. Our results show that, in the context of 
nucleic acid hairpins, precise measurement of the sequence-dependent folding landscape of 
RNA is sufficient to obtain good estimates of the /-dependent hopping rates in the absence of 
handles. It suffices to measure F eq (z sys ) ~ F eq (z m ) ps F° q (z m ) at / = f m using stiff handles, 
while F eq (z m \f) at other values for / can be obtained by tilting F eg (z m \f m ). The accurate 
computation of the hopping rates using F eq (z m ) show that z m is an excellent reaction coordinate 
for nucleic acid hairpins under tension. Further theoretical and experimental work is needed to 
test if the proposed framework can be used to predict the force dependent hopping rates for 
other RNA molecules that fold and unfold through populated intermediates. 



METHODS 



RNA hairpin: The Hamiltonian for the RNA hairpin with N nucleotides, which is modeled 
using the self-organized polymer (SOP) model 35j, is 



H 



kRl 



SOP 



N-l 

5> 



i=l 



1 



Rl 



N-3 N 
i=l j=i+3 



i -J 
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N-2 



'.J 



'.J 



12 



r o \ 6 

21 



i.j 



A 



i.j 



i=l 



a 



i,i+2 



(2) 



where r^j = |rj — rj\ and r°j is the distance between monomers i and j in the native structure. 
The first term enforces backbone chain connectivity using the finite extensible nonlinear elastic 
(FENE) potential, with k w 1.4 x 10 4 pN-nm" 1 and R = 0.2 nm. The Lennard- Jones interaction 
(second term in equation (j2J)) describes interactions only between native contacts (defined as 
T°j < 1.4 nm for \i — j\ > 2), with Ajj = 1 if monomers i and j are within 1.4 nm in the 
native state, and Ajj = otherwise. Non-native interactions are treated as purely repulsive 
(the third term in equation (j2J)) with a = 0.7 nm. We take = 4.9 pN-nm and e\ = 7.0 pN-nm 
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for the strength of interactions. In the fourth term, the repulsion between the i th and (i + 2) th 
interaction sites along the backbone has a* = 0.35 nm to prevent disruption of the native helical 
structure. 

Handle polymers: The handles are modeled using the Hamiltonian 

, N-l N-2 

Hhandles = ~ ~~ r o) 2 ~~ ' *i+l,i+2- (3) 

i=l i=l 

The neighboring interaction sites, with an equilibrium distance r = 0.5 nm, are harmonically 
constrained with a spring constant ks ~ 1.4 x 10 4 pN-nm -1 . In the second term of eq. ([3]), the 
strength of the bending potential, k^, determines the handle flexibility. We choose two values, 
k\ =7.0 pN-nm and k^ =561 pN-nm to model flexible and semiflexible handles respectively, and 
assign k& = 35 pN-nm to the junction connecting two ends of the RNA and the handles. We 
determine the corresponding persistence length for the two k& values as l p = 0.6 and 70 nm (see 
SI text). The contour length of each handle is varied from N = 5 — 200. 
Generalized Rouse model (GRM): The Hamiltonian for the GRM (Fig. 3a) is 

o N h 2N h +N 

i=l i=N h +N +l 
g N h +N 

+ ^2 Yl (^+i-n) 2 + PV[r N - Nh+1 -r Nh+l ], (4) 

i=N h +l 



where 




V[v] = ' " . (5) 



The first two terms in equation (j3J) are the discrete connectivity potentials for the two handles, 
each with Nh bonds (Nh + 1 monomers), and with Kuhn length b. The mechanical force f in the 
third term is applied along the z direction, with |f | = f m = 15.4 pN. We also fix the 5' end of the 
system with a harmonic bond of strength k = 2.5 x 10 4 pN-nm -1 in the fourth term of eq. (j3J). 
The fifth term mimics the RNA hairpin with iVo bonds and spacing a = 0.5 nm. Interactions 
between the two ends of the RNA hairpin are modeled as harmonic bond with strength k pa 0.54 
pN-nm -1 that is cut off at c = 4 nm (eq. (jSJ)). When |R m | exceeds 4nm, the bond is broken, 
mimicking the unfolded state. 
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The free energies as a function of both R m « z m and R sys ~ z S3/s are most easily determined 
in the continuum limit of the Hamiltonian in equation (j3J), with Ylf=i ~* Jo ds. Because of 
the relatively large value of the external tension (f m ^> ksT/lp), we can neglect transverse 
fluctuations without significantly altering the equilibrium or kinetic properties of the GRM 
The refolding time, t™ of the RNA mimic (Fig. 3a), which is the WF closure time 
can be determined by numerically diagonalizing the Rouse-like matrix with elements 
Mij = l^H/STiSTj. 
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FIG. 1: A schematic diagram of the optical tweezers setup used to measure the hairpin's folding 
landscape, a. Two RNA/DNA hybrid linkers are attached to the 5' and 3' ends of the RNA hairpin, 
and a constant force is applied to one end through the bead. b. Ensemble of sampled conformations 
of the H-RNA-H system during the hopping transitions obtained using L=25 nm and l p = 70 nm. 
The illustration is created using the simulated structures collected every 0.5 ms. An example of the 
time trace of each component of the system, at / = 15.4 pN, is given L for both linkers is 25 nm. 
z m {= Zy — zy) and z sys (= z p — z Q ) measure the extension dynamics of the RNA hairpin and of the 
entire system respectively. The time averaged value zxit) = j Jq dTz sys (T) for the time trace of z sys is 
shown as the bold line. The histograms of the extension are shown on top of each column. 
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FIG. 2: The free energy profiles, F eq (z sys ) (dashed line in blue) and F eq (z m ) (solid line in red), calculated 
using the histog rams obtained from the time traces z S y S {i^j and z m (t^) for varying L and l p . F eq (^z S y S ^ 
and F eq (z m ) for a given l p and different L are plotted in the same graph to highlight the differences. 
The intrinsic free energy F° q {z m ), the free energy profile in the absence of handle, is shown in black. 
The condition that produces the least deviation {l p = 70 nm, L = 5 nm) and the condition of maximal 
difference (l p = 0.6 nm, L = 25 nm) between F eq (z m ) and F eq (z sys ) are enclosed in the magenta and 
orange boxes, respectively. 
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FIG. 3: Free energy profiles for the GRM. a. A schematic diagram of the GRM, showing the number 
of monomers (Nq and N^) and Kuhn lengths (a and b) in each region of the chain, the harmonic 
interaction between the ends of the RNA mimic, and the external tension, b. The free energy profile 
for a fixed b(= a/3) and increasing Nh as a function of R sys ~ z sys . The barrier heights decrease and 
the well depths increase for increasing N^. c: The free energy profile for fixed = 20 and varying b. 
The barrier heights decrease and the well depths increase for increasing b. In both b and c, the profiles 
are shifted so that the positions of the local maxima and minima coincide with those of the intrinsic 
free energy (with Nh = 0). 
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FIG. 4: a. The free energy profile for P5GA with L = nm. b. The transition times at / = f m , 
obtained using z m (t) (filled symbols) and z sys (t) (empty symbols). The ratio T u{F)l T u{F) I, which 
shows that z sys (t) mirrors the hopping of P5GA. c. Folding rate k™ (L) / k™ (0) as a function of L for 
varying b, using the GRM. The plots show b/a =1, 1/2, 1/5, and 1/10. The simulation results for 
P5GA are also shown as symbols, to emphasize that the GRM accounts for the hopping kinetics in the 
H-RNA-H system accurately. 
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FIG. 5: a. Comparison of the measured free energy profiles (symbols) with the shifted free energy 
profiles f3F° q {z m \f m ) — f3(f — f m ) • z rn . b. Folding and unfolding times as a function of force / = 14 
pN < f m , f = 15.4 pN « f m , and / = 16.8 pN > f m . is obtained from the time trace in 



Figure 2B in Ref [35j at each force, while T^im is computed using the tilted profile l3F° q (z m \f) 
f3F° q (z m \f m ) — /?(/ — fm) ■ z m in equation ([T]). c. Folding and unfolding times using the GRM. Symbols 
are a direct simulation of the GRM (error bars are standard deviation of the mean) . The solid lines are 
obtained using the Kramers theory (equation {!])). We choose Djj ~ 3Dq, so that that the simulated 
and Kramers times agree at / = f m . The position of each basin of attraction as a function of force for 
the GRM is given by z v w N a 2 (3f/3 and z F w N a 2 (3f/{3 + 2N a 2 (3k). 
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SUPPORTING INFORMATION (SI) 



Fluctuations of the handle polymer under tension : Using the force clamp simulations 
of the RNA hairpin in the presence of handles, the dynamics of the fluctuations of the handle can 
be independently extracted by probing the time-dependent changes in the 5' and 3' ends of the 
RNA molecule. The distribution of the longitudinal fluctuations (z^ = z^ — z and = z p — zy) 
and the dispersion in the transverse fluctuations (x^ or y 5 > and xy or yy) are shown in SI Fig. 
6. Assuming that the force f = f\\e\\ + f±e± is decomposed into / = |f | ~ f\\ + f]_/2f\\ + 0(fj_) 
(with /|| ^> and using the partition function of the worm-like chain polymer under tension, 
Z = J d{le~( HwLC / kBT ~f'^ kBT ^ jlj], one can express the longitudinal fluctuation as 

, 9, d{R\\) \ LL f or / < MI 

(SRh = k B T^-j^ ~ I p , , 6 

K " d f [Ll^if/knT)-^ for/»*f 

where the force extension relations of a worm-like chain R\\/L m fl p /k B T for / < k B T/l p and 
R\\/L « 1 - ^k B T/Al p f for / > k B T/l p are used for small and large forces, respectively [2|. 
These results are consistent with the fluctuations observed in the simulations. When / < k B T/l p , 
the transverse fluctuations are independent of the force, and are determined solely by the nature 
of the linker. For / = J m ~ 15.4 pN, the tension is in the regime that satisfies / > k B T/l p for 
both values of l p used in the simulations, and the longitudinal fluctuations (5R?i) decrease as 
the stiffness of the polymer increases for all L (SI Fig. 6a). The distribution of the extensions 
coincide for both the 3' and 5' ends of the handles for all L and l p (SI Fig. 6a). This suggests that 
the constant force applied at the point z p propagates uniformly throughout the whole system. 
The transverse fluctuations are given by 

(SR ± ) - (k B T) -^-1,,=,,^ - — <*„> ~ , _ ™ for (7) 



/ V 2 V l p f . 

The transverse fluctuations also decrease as /(> k B T/l p ) increases, with a different power. It 
is worth noting that the transverse fluctuations, which also increase as L increases, are nearly 
independent of the handle stiffness if / > k B T/l p . The standard deviations of distributions are 
plotted with respect to the contour length at each bending rigidity (SI Fig. 6b). The fit shows 
that a ~ 0.28 x L ' nm and a ~ 0.30 x L 1 ' 2 nm for l p = 0.6 nm and l p = 70 nm, respectively, 
which is consistent with the analysis in equation [71 
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Determination of the persistence length of the handles : In order to determine the 
persistence length of the handles, we numerically generated the end-to-end distribution function 
P(R) of the free handles in the absence of tension, and fit the simulated distribution to the 



analytical result 



, . AirCp 2 ( a \ . . 

Pwlc{R) = L(1 _ p2)9/2 e*P [~Y^ ) ^ 

with p = R/L and a = 3L/4l p . The normalization constant 

C =[vr 3 / 2 e- Q a 3 / 2 (l + 3a" 1 + Iha^/A)]" 1 ensures ft dR P W lc{R) = 1- 

Quantifying the synchronization of z sys (t) and z m (t) : The origin of the small discrep- 
ancy between F eq (z sys ) and F eq (z m ) when the handles are flexible can be found by comparing 
z S ys{t) with z m (t) for the two extreme cases in Fig. 2 of the main text. To quantitatively express 
the synchronization between z sys (t) and z m (t), we defined a correlation function at each time t 
using 

n(A z sys(t) - Z syg Z m (t) - Z^f 

C W jts x Its W 

Sys 771 

where z^ y s s and z^ s are the positions of the transition states determined from F eg (z sys ) and 
Feq(z m ) respectively. If C(t) > at time t, both z sys (t) and z m (t) are in the same basins 
of attraction, i.e. the status of z m (t) is correctly detected by the measurement through the 
handles. If C(t) < 0, then the information of z m (t) is lost due to fluctuations or the slow 
response of the handles. The near perfect synchronization between z sys (t) and z m (t) for l p = 60 
nm and L = 5 nm are reflected in C(t) > for almost all t. Thus, when the handles are stiff, 
z m (t) ~ Zgy S (t) — 2L, which implies that z sys {t) faithfully reflects the dynamics (z m (t)) of the 
hairpin. In contrast, with l p = 0.6 nm and L = 25 nm, z m (t) can not be determined from 
z S y S (t) using z sys (t) — 2 x L ^ z m (t). The amplitudes of z sys (t) are typically larger than that 
of z m (t), leading to C(t) < occasionally (shown by an arrow on the right plot in SI Fig. 7). 
The histograms of P{C) for the two extreme cases show that the dynamics between z sys (t) 
and z m (t) are more synchronous for the rigid and short handles (0.0 < P{C) < 0.5) than for 
the flexible and longer handles (—0.05 < P{C) < 0.2) (see the graph at the bottom in the SI 
Fig. 7). The finding that short and stiff (L/l p ~ 0(1)) handles minimize the differences between 
F e q(z m ) and F eq (z sys ) is related to the tension-dependent fluctuations in the linkers. 
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FIG. 6: Fluctuations of the handles with varying lengths and flexibilities at / = 15.4 pN. a. Longi- 
tudinal fluctuations of the handle attached at the 5' and 3' sides of the RNA hairpin, b. Transverse 
fluctuations are fit to a Gaussian distribution, and the standard deviation (a) is plotted as a function 
of the contour length and flexibility. 
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FIG. 7: The time traces, z sys (t) and z m (t), for two extreme cases that produce the free energy profiles 
in the magenta and orange boxes in Fig. 2 of the main text, are overlapped to show the differences. For 
l p = 70 nm and L = 5 nm, both the phase and amplitude between z sys (t) and z m (t) coincide throughout 
the time series, while for l p = 0.6 nm and L = 25 nm the amplitude of z sys (t) are larger than z m (t) and 
the phase between z sys (t) and z m (t) is occasionally offset from one another. The correlation measure 
C(t) quantifies the synchrony between z sys {t) and z m (t) at time t. The histograms of C(t) show that 
the time trace for l p =70 nm and L = 5 nm is more synchronized than the one for l p = 0.6 nm and 
L = 25 nm. 
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FIG. 8: The survival probabilities Pp(t) and Pu(t) are fit to a single exponential function to calculate 
tjj and tf- For L=0 nm, iy = 2.9 ms and Tp = 1.9 ms. For L=20 nm, tjj = 5.0 ms and Tp = 12.1 
ms. The quality of the fits for L = 20 nm (dashed lines) is not as good as for L = nm. The survival 
probabilities show lag phases for both unfolding and refolding. 
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FIG. 9: Time traces of molecular extension under tension / =14.0, 15.4, and 16.8 pN, and corresponding 
distribution P(z m \f) at each force. The distributions are converted to the free energy profile in Fig. 5a 
by using F° {z m ) /ksT = — \ogP(z m ). Note that the hairpin is pinned in the UBA at / = 16.8pN 
(> fm) with infrequent transitions to the NBA. Just as in experiments accurate measurement of 
Feq(zm) 1S possible only at / ~ f m , where multiple hopping events between the NBA and UBA can be 
observed. 



26 



